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ABSTRACT 



Aims. The aim of this work is a detailed analysis of transit light curves from TrES-1 and TrES-2, obtained over a period of three to 
four years, in order to search for variabilities in observed mid-transit times and to set limits for the presence of additional third bodies. 
Methods. Using the lAC 80cm telescope, we observed transits of TrES-1 and TrES-2 over several years. Based on these new data and 
previously published work, we studied the observed light curves and searched for variations in the difference between observed and 
calculated (based on a fixed ephemeris) transit times. To model possible transit timing variations, we used polynomials of different 
orders, simulated 0-C diagrams corresponding to a perturbing third mass and sinusoidal fits. For each model we calculated the 
residuals and the False Alarm Probability (FAP). 

Results. For TrES- 1 we can exclude planetary companions (> 1 M^) in the 3:2 and 2: 1 MMRs having high FAPs based on our transit 
observations from ground. Additionally, the presence of a light time effect caused by e. g. a 0.09 Mq mass star at a distance of 7.8 
AU is possible. As for TrES-2, we found a better ephemeris of T, = 2 453 957.635 12(28) + 2.4706101(18)xEpoch and a good fit for 
a sine function with a period of 0.2 days, compatible with a moon around TrES-2 and an amplitude of 57 s, but it was not a uniquely 
low x~ value that would indicate a clear signal. In both cases, TrES-1 and TrES-2, we were able to put upper limits on the presence 
of additional perturbers masses. We also conclude that any sinusoidal variations that might be indicative of exomoons need to be 
confirmed with higher statistical significance by further observations, noting that TrES-2 is in the field-of-view of the Kepler Space 
Telescope. 

Key words, planetary systems - Methods: N-body simulations - Techniques: photometric 



1. Introduction 

In 1992, the first exoplanets around the pulsar PSR B1257+12 
were discovered by detect ing anomalies in the pulsation period 
dWoIszczan & FraiH Il992h . Similarly, by studying variations 
in the time of transit occurrence, transiting exoplanets give 
another possibility to find additional companions, even down 
to Earth masses. There are different mechanisms causing these 
variations. For one the gravitational influence of a perturb- 
ing body can alter the or b ital period of the transiting planet 
dHoIman & Murravl 120051: lAgoI et al.1 \2O05). This effect is 
particularly sensitive to detect additional bodies in mean-motion 
resonances with the transiting planet or to detect moons around 
that planet. For another, a perturbing mass in an orbit larger 
than the transiting planet can cause the "star - transiting planet" 
system to wobble around the barycenter and alter the observed 
period i city, something that is known a s the light-time effect 
dlrwiiil \195^ ISchneider & Dovia [1991 iDovIe & Deeal l2004t 
ISchneideT 2005 ). These anomalies are reflected in the times 
of the transit occuiTence. Hence, using a linear ephemeris 
and several observations of transits, it is possible to get the 
difference 'O-C between the observed transit times and the 
calculated ones from the ephemeris. From this difference we 
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can learn about perturbations on strict orbital periodicity due to 
a possible companion object which has not been detected yet. 
Motivated by the possible detection of low-mass companions 
around transiting planets, we started a lon g-term c ampaign to 
observe tran sits of the exoplanets Tr ES-1 jAlonso et al. 200^ 
and TrES-2 dO'Donovan et alj|2006l) over several years. Here 
we present the results of our observations and a study of the 
transit timing variations that were found. 

Searches for planet-mass objects from the light time effe ct were 
first p roposed for low-mass eclipsing binaries by Do vle et al.l 
(Il998l) and have been performed in-depth o n the system CM- 
Draconis, with ambiguou s results to date (iDeeg et al.l l2000l 
120081: [Morales et alJl2009h . They were soon followed by several 
studies of timing effects in transiting planet systems: TrES-1 
(Steffen & Agol 2005), HD 189733 (Mifler-Ricci et al. 2008^), 
HD 209458 (Agol & Steffen 2007; Miller-Ricci et alj (20081, 
GJ 436 tAlonso et al., ,2008; ,Bean & Seifahit ,200i) and for 
CoRoT-lblBeanI 12003)" In all cases the authors could only 
constrain the parameters of a potential second planet. 
The first search for variability in the t ransit times of TrES-1 
was conducted by lSteffen & Agoll (l2005l) . They fitted perturbing 
planet s in different orbi ts , usin g transit timing o bserva tions 
from ICharbonneau et al.1 (l2005h . ISteffen & Agol (|2005[) did 
not fi nd any strong evidence of a third object. IWinn et aP 
(l2007i) observed three additional transits of TrES-1; they found 
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that the transit times of these three observations occurred 
progressively later. However, as mentioned by the authors, 
these three measureme nts were no t enou gh to make any firm 
conclusions. Recently, iRaetz et aP (l2009b ) reported a transit 
timing analysis of TrES-1 with the only result of being an im- 
proved ephemeris. Previo us studies of tr a nsit tim ing s for TrES-2 
have be en performed by 'Holm an et al.l (l2007h and iRaetz et aP 
(l2009ah . neither finding e vidence fo r an unseen perturbing 
planet. Very recently Misl is & Schm itt (2009) analysed transit 
time durations between 2006 and 2008, indicating a possible 
change of inclination. 

In Section|2]we describe the procedure to fit our observations to 
a transit model and to obtain the observed mid-transit times and 
their error estimates. These values are then used in Section[3]to 
interpret the O-C diagram by means of different models. Finally 
we discuss the results in Section H] 



2. Light curve fitting 

All observations reported here were performed with the lAC 
80cm telescope (IAC-80) at the Teide Observatory, Tenerife. We 
used the same observing and analysis procedure for both TrES- 
1 and TrES-2. For further investigation described here, we only 
used observations where a complete transit light curve was ac- 
quired. A detailed description of the analysis leadi ng from the 
telesco pe images to the light curves can be found in lRabus et al.l 
( 1200 9V where a study of flux variations during transits due to 
starspots or additional transiting planets for TrES-1 was pre- 
sented. We decorrelated the light curves against the airmass by 
subtracting a quadratic polynomial fit. Due to a slight defocus- 
ing of the telescope during observations, the centroids are not 
well defined and de-correlation against the target's detector po- 
sition is not applicable. However, the de-focusing and spread of 
the flux over several pixels will lower the noise correlated with 
the target's detector position. 

For TrES-1 we obtained eight useful transit observations (Table 
|2]l over a period of three years and for TrES-2 five observations 
(TableO over a period of two years. In order to measure the ob- 
served mid-transit times, we first created a template of the transit 
event from a folding in phase and by binning (6-point bins, with 
a mean size of 66 s for TrES-1 and 90 s for TrES-2) of all obser- 
vations of the respective transiting planet (Figure [1). The stan- 
dard deviations inside the individual bins were ~ 1.63 mmag 
for TrES-1 and ~ 1.79 mmag for TrES-2, whereas the stan- 
dard deviation outside the transit part of the binned light curve 
is 0.7 mmag for TrES-1 and 0.8 mmag for TrES-2. 
We used these phased and binned light curves to creat e a model 
of the transit light curve, using the formalism from GimenezI 
(|2006|) and the simplex-downhill fitting algorithm dPress et alj 
Il992h . We fitted for the planetary and stellar radii ratio, k, the 



sum of the projected radii, rr, and orbital inclination, /, while 
we kept the eccentricity fixed at zero and th e limb darkening co - 
efficient fixed, as obtained by the tables of IClaret et al.l(ll995h . 
assuming quadratic coefficients. The best-fit models are plotted 
over phased data in Fig. [1] and the best-fit parameters are shown 
in Table[T] We can see that our light curve model param eters are 
consistent with the parameters from lSouthworthI (|2008|) for both 
TrES-1 and TrES-2. 

We then shifted in time our model with the fixed best-fit param- 
eters against each individual observed light curve and calculated 
the;t'^ residuals for each shift. The fit with the minimum value 
gave us the observed mid-transit time. We estimated the timing 
error within a 68 % (1-cr) confidence interval of the values. 



given by the range where the;t'^ residuals increases by 1 over the 
minimum value. We also calculated the timing precision of each 
individual light curve, 5,, by propagatin g their respective pho to- 
metric precision, 8l, using the equation dPovle & Deeg||2004b : 



L(f,_i)-L(f,+i) 
2Af 



(1) 



where L(f,) is the stellar brightness at f,. At is the cadence and 
the sum goes over all brightness values within the eclipse event. 
Comparing the estimated error, based on the 1-cr interval, with 
the error calculated using Equation [T] we found a good agree- 
ment; both error values are indicated in Tables|2]and[3] Generally 
the estimated 1-cr errors are higher than the propagated errors, 
but in four cases the error based on the 68 % confidence interval 
were lower. Our timing errors are possibly underestimated due 
to correlated noise. However, we do not modify them in order 
to stay consistent with the timing measurements and associated 
errors that we took from the literature. 

Finally we subtracted the observed mid-transit times from 



Table 1. Compa rison b etween best-fit values of this work and 
parameters from lSouthw orth (2008) for transit light curve mod- 
els of TrES-1 (upper values) and TrES-2 (lower values). 



Parameter 


This work 


Southworth (2008) 


TrES-1 






k 


0.1350 


0.1381 ±0.0014 


rr 


0.1104 


0.1097 ±0.0022 




88.67 


88.67 ± 0.71 




TrES-2 






k 


0.1260 


0.1268 ±0.0032 


rr 


0.1462 


0.1460 ±0.0042 


in 


83.70 


83.71 ±0.42 



the calculated ones, obtaining the O-C values. The calculated 
mid-transit times for TrES-1 were obtained from the ephemeris 



Table 2. Overview of TrES-1 transits used in this work. 



Epoch 


observed 


O-C 


l-o- 


Calculated 


Source 




transit time 


[s] 


error 


timing 






HJD-2450000 




[s] 


precision [s] 




-4 


3174.6864 


60 


35 




1 


-1 


3183.7752 


-63 


43 




1 





3186.8061 


9 


26 




1 


20 


3247.4075 


2 


35 




1 


124 


3562.5352 


14 


20 


13 


2 


126 


3568.5952 


-6 


22 


25 


2 


234 


3895.8430 


-24 


16 




3 


235 


3898.8734 


8 


12 




3 


236 


3901.9037 


28 


16 




3 


254 


3956.4445 


-15 


14 


20 


2 


255 


3959.4760 


107 


17 


23 


2 


357 


4268.5418 


-44 


14 


17 


2 


358 


4271.5716 


-66 


17 


16 


2 


359 


4274.6028 


32 


18 


17 


2 


386 


4356.4142 


-24 


9 




4 


387 


4359.4443 


-18 


13 




4 


388 


4362.4742 


-31 


17 




4 


490 


4671.5412 


-56 


26 


29 


2 



Sources: 1) Charbonneau et al. (2005), 2) This work, 3) Winn et al. 
(2007). 4) hrudkova et al. (2009) 
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Table 3. Same as Table |2l for TrES-2. 



Epoch 


observed 


O-C 


1-0- 


Calculated 


Source 




transit time 


[s] 


error 


timing 






HJD-2450000 




[s] 


precision [s] 




4 


3967.5180 


92 


37 


22 


1 


13 


3989.7529 





25 




2 


15 


3994.6939 


-15 


27 




2 


34 


4041.6358 


-10 


26 




2 


140 


4303.5209 


-72 


26 


20 


1 


142 


4308.4613 


-169 


39 


22 


1 


274 


4634.5828 


-184 


26 


24 


1 


276 


4639.5232 


-257 


27 


21 


1 



Sources: 1) this work, D Holman et al. (2007) 



Tc = 2453 186.8060 + 3.0300737x Epoch dWinn etal.ll2007h 
and for TrES-2 we us ed T,. = 2453 957. 6348 +2.47062 Ix Epoch 
jHolman et alj|2007l) . 

In the following transit timing analy sis we included seve r al pub - 
lished m id-transit observations, by Charbonn eau et all (l2005h . 
IWinn et al. (2007) and Hrudkova et al. (2009) for TrES-1 (Table 
m, and for TrES-2 by Holman et alj (12007 ,) (Tabled. We con- 
sidered only O-C times with errors below 60 sec, which led 
to the rejection of some O-C values from ICharbonneau et all 
( I2005h . In the O-C residuals of TrES-1 we also removed two 
outliers at Epoch 2 55 and 358 which have been identified by 
iRabus et"an (l2009l) as transits with possible starspots. Several 
further transit times of TrES-2 and TrES-1 were recently re- 



ported by 'Raetz et al.' ('2009a"^, respectively, and for TrES-2 
by Mishs & Schmitt (2009) , but we did not include these in our 
study, since they had individual errors and an internal scatter sev- 
eral times larger than the data included in this study. The tran- 
sit mid-time corresponding to Epoch of TrES-2, reported by 
lO'Donovan et alj ([2006) has also been removed due to its high 
error of >60 s. We obtained for TrES-1 an O-C diagram spanning 
four years, with 16 points (Figure |2| and for TrES-2 a diagram 
spanning two years, with 8 points (Figure[3]l. 



3. Transit timing analysis 

Two mechanisms have been brought forward that may cause 
deviations of transit times from strict periodicity. For one the 
light time effect and for another the influence of a third body on 
the transiting planet's intrinsic periodicity. 



3. 1 . Search for the light time effect 

We first consider the light time effect, where the barycenterof the 
binary system, in our case the system "star - transiting planet", is 
offset against the barycenter formed by an additional third body. 
This will cause the light to travel a longer or shorter distance to 
the observer and hence the observer will see a different period, 
where the intrinsic period is unknown. 

For the analysis of the O-C times, we fitted three polynomials 
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Epoch 

Fig. 2. Difference between calculated (based on ephemeris from 
IWinn et all (|2007|) ) and observed mid-transit times for TrES-1. 
Filled dots are O-C values obtained w ith the IAC-80, squares are 
taken f rom lCharbonneau et a L ( 2005) , triangle s from Wi nn et aTl 
(12001) and asterisk are from Hrudkova et al.l (12009). The lines 
show different polynomial fits, where the solid black line indi- 
cates a linear fit, the long-dashed (red) line a quadratic and the 
short-dashed (blue) line a cubic polynomial. The dotted (red) 
lines show the fits corresponding to the variation of the quadratic 
term within l-cr confidence limits. 



200 p 
100 ^ 




50 100 150 200 250 300 

Epoch 

Fig. 3. Simil a r as Figure |2] (based on ephemeris from 
iHolman et al.l (l2007h ) including only polynomials fits and cor- 
responding l-cr confidence limits, but for TrES-2. Filled dots 
are O-C values obtained w ith the IAC-80 and squares are taken 
from lHolman et alj (l2007h . 



of different orders to the O-C values, namely a linear, OC fi, — 
kq + Eku a parabolic, OC//, - a-q + Eki + E^K2, and a cubic 
polynomial, OC/i, = kq + Ek\ + E^K2 + E^Ky,, where E is the 
Epoch number and k the fitted polynomial coefficients. We also 
analyzed the case of fixed original ephemeris (OC/i, = 0) and 
the case of maintaining the original period and fitting only for a 
constant offset in O-C (OC/,, = ko). Table|4]shows the obtained 
best fit coefficient for TrES-1 and TrES-2 respectively, and Figs. 
|2] and [3] show a plot of the respective polynomials. In order to 
quantify the improvement of the different fits, we calculated the 
False Alarm Probability (FAP), which is the significance level 
of the fit quality improvement and indicates the probability of 



making a type I error. Therefore, we first calculated the F-values 
according to equation: 

^_ (yf-^^)/(vi-v2) 

where are residuals from the lower order fit, X2 residu- 
als from the higher order fit, and V] 2 are the corresponding de- 
grees of freedom (Table |4|i. We then used the F-value to calcu- 
late the significance using the IDL-routine MPFTEST from the 
Markwardt IDL library, which gives the probability for a value 
drawn from the F-distribution to equal or exceed the given F- 
value. For TrES-1, FAFs were calculated against the original 
ephemeris, against the offset-only case and against the linear 
polynomial. We see that the original ephemeris is unlikely to be 
the best solution with the linear and quadratic fits being the most 
likely descriptions, having low residuals and a low FAP. For 
TrES-2 a linear trend in the O-C residuals (Fig. O is apparent. 
The statistical analysis shows no clear preference between a lin- 
ear or a quadratic polynomial (Table |4]i. The distribution of the 
observations of TrES-2 into three groups, acting as pivot points, 
does not support fits of orders higher than the quadratic one. 
Therefore, we fitted no cubic polynomial to the TrES-2 data. 
As for the linear coefficients /cq and ki these are without relevant 
physical meaning, but indicate a slightly different ephemeris 
than used, where a-q is an offset of the mid-transit time at epoch 
and Ki a correction to the period of the ephemeris. However, 
we can give a physical meaning to the quadratic coefficient. The 
quadratic term {K2), gives the s ystem's acceleration along the line 
of sight by using the equation (iDeeg et al...2008i) : 



where c is the speed of light and P is the observed period. 
In order to estimate the error of the quadratic term, we stepped 
through different values of K2, fitting for kq and ki . The K2 values 
which increased the best-fit residuals by 2.3 gave the l-cr con- 
fidence limits (see Table|4]i; dotted lines in Figs.|2]and[3]show the 
corresponding fits at these limits. Our best fit parabola model for 
TrES-1 gave a quadratic term of a:2 = -3.6±3.4x lO"'' s. Solving 
Eq.[3]we obtained an acceleration of ay = -3.2 + 3.0x 10"^ m s"- 
and — = -1.1 + 1.0 X 10"'^ s"'. The cubic term /C3 would indi- 
cate a constant change in acceleration; however the cubic fit is 
also less likely than the quadratic one and will not be further dis- 
cussed. 

For the quadratic solution of TrES-2 we obtained a best-fit value 
of K2 smaller than its error. This high error is consistent with 
a low significance of the quadratic solution, therefore we sup- 
port the linear case, with a l-cr upper limit for accelerations of 

fl|l < 3.2 X 10"^ ms--. 

Considering the clear linear trend of O-C times in Figs. |2]and[3] 
and using the coefficients from the linear fit, we indicate here an 
improved ephemeris for TrES-1, given by: 

=2453 186.8061 1(16) -f 3.0300728(6) X Epoch, (4) 

where the values in parenthesis give the uncertainty in the last 
two digits and for TrES-2: 



=2453 957.63512(28) -H 2.4706101(18) X Epoch. (5) 
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Fig. 8. Simulated O-C diagram of TrES-2 and the best fit perturber with a mass of 18 Me and an orbit of 0.051 AU. The red dots 
coiTespond to the observations and the squares to the coiTesponding simulated O-C values. 



3.2. Search for a perturbation of the intrinsic period 

The other cause for transit timing variations could be a per- 
turbation of the intrinsic planet period P' due to the grav- 
itational influence of a third body on the transiting planet. 
Regarding such perturbations, there exist no analytical equa- 
tions that describe the gravitational influence on a transit- 
ing planet due to an undetected third body. Generally, N- 
body simulations are used, iterating over a large orbital pa- 
rameter space for a maximum possible mass rang e, see e. g. 
[Holman & Murrav (2005); Steffen & Agol (2005): lAgol et al.1 
(2005); Agol & Steffen (2007); iMiUer-Ricci et al] (l2008allbl) ; 
iNesvorny & Morbidelli ( 2008) . 

In order to find the parameter space of third bodies compati- 
ble with the observed transit times, we created a numerical 2- 
dimensional simulation of a three body system by integrating 
over the equations of motion. We considered the problem in 
two dimensions, assuming that the orbits of the exoplanetary 
system and its respective perturber are co-aligned. We further 
considered the problem in a helio-centric frame, with the star 
in the center of the coordinate system, meaning that we neglect 
the light-time effect in this context. To integrat e the equation of 
motion we used the Burlisch-Stoer algorithm (.Press et alJI 19921) 
with a 1 second time step and an accuracy of 10"'". 
We considered the problem for an inner and an outer per- 
turber separately and neglected the 1 : 1 mean motion resonance 
(MMR). Each simulation started with a perturber 0.005 AU away 
from the transiting planet, with an initial zero eccentricity and 



with a phase shift towards the transiting planets between 0° and 
315° in steps of 45°. We advanced the simulations in orbital steps 
of 0.001 AU and used 100 mass steps from 1 to 100 M® for 
the perturbing object. We simulated 1000 transits, correspond- 
ing to approximately eight years coverage for TrES-1 and seven 
years coverage for TrES-2. We used a mass o f 0.61 M j for TrES - 
1 and a mass of 0.88 Mq for its host star dTorres et al.ll2008l) . 
and a mass of 1 .2 M / for TrES-2 and 0.98 Mq for its host star 
(ISozzetti et al.l 12007 *). After obtaining the simulated O-C dia- 
grams, we applied a Fourier transformation to the synthetic O-C 
diagrams (iMiller-Ricci et al.t2 008a b) and derived the maximum 
obtained O-C amplitudes as a function of the perturber's semi- 
major axis; showing them for masses of 10, 20 and 40 in Fig. 
m where also the MMRs are indicated by vertical lines. Figs. |5] 
and|6]show some examples of synthetic O-C diagrams for TrES- 
1. 

To establish potential third body orbital parameters, we fitted 
for each perturber distance the observed O-C diagram to all syn- 
thetic ones that had been generated for masses of 1 - 100 
(Fig. |7]i. We left three parameters open in the fits: one param- 
eter to shift the O-C diagrams in Epoch, coiTesponding to find 
the best moment for the first real observation (E=0) within the 
simulated data, and two parameters describing linear deviations 
between simulated and modelled O-C diagrams, which have no 
relevant physical meaning. 

For the best-fit masses, we see a general trend to higher masses 
for perturbers closer to the host star and at larger semi-major 
axis, and we can also identify the mean-motion resonances. 






0.0100 0.0145 0.0190 0.0235 0.0280 
semi-major axis [AU] 



10"' I I 

0.044 0.058 0.072 0.086 0.100 
semi-major axis [AU] 



Fig. 4. Maximum O-C amplitudes of a 2 dimensional N-body 
simulation for the system TrES-1 (upper two panels) and TrES-2 
(lower two panels) due to a perturber during 1000 transits plotted 
against the semi-major axes of the perturber The colored solid 
lines correspond to different perturbing masses, namely 10 
(red), 20 M© (blue) and 40 M© (green). Left side: inner perturber, 
right side; outer perturber. Vertical lines indicate the MMRs. 



where the best fit for a consistent perturber indicates lower 
masses (vertical lines in Fig. At larger orbits for the per- 
turber, the best-fit perturber's mass increases and it increases 
more steeply for TrES-2 than for TrES-1. This is due to TrES-2 
having a higher mass than TrES-1 and being closer to the cen- 
tral star. The high peaks at the 2:3 and 2: 1 MMRs for TrES-1 
reveal that a >1 M© planet in that configuration might have been 
detected in our data set. 

Again, in order to quantify the improvement of the simulated O- 
C fits in the lower MMRs and the best-fit simulated O-C against 
the linear polynomial, we calculate the FAPs. In these cases it is 
better not to use Eq.|2]to calculate the F-value, since this equa- 
tion is better suited for cases where one expects small changes 
for additional parameters, whereas for the simulated O-C values, 

we expect big changes. Therefore, we used F - ^ly^^ to calcu- 
late the F-value for the MPFTEST-routine. In Table |4l we show 
the and FAPs for some special cases of our simulated O-C 
values. 

For TrES-1 we found that most of the;^'^ values for the simulated 



Fig. 5. Example synthetic O-C diagrams for TrES-1 and a per- 
turber of 30 Me at 0.053 AU and with different phase shifts be- 
tween them. Note that the appearance of several lines is due to 
an aliasing effect between consecutive transits. For a given tran- 
sit number there is only one O-C time. 
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Fig. 6. Same as Fig. |5] but for a perturber distance of 0.062 AU 
(2:1 MMR). 



O-C models were below the one of the linear fit, but based on 
the available data, we can exclude the case of 2:3 and 2: 1 MMR 
for planets more massive than 1 M© with high FAPs. We found 
a best-fit for a perturber with a mass of 2 M© and a semi-major 
axis of 0.05 AU, but this peak is not outstandingly low. Similarly, 
for TrES-2, all values for the O-C models are below the one 
from the linear fit. We can identify some low peaks, with the 
lowest at 0.049 AU (x^ - 0.3) for a perturbing mass of 18 M^. 
Fig. [8] shows for TrES-2 the O-C simulation for the lowest x^ 
with the observed O-C residuals over-plotted. Given that none 
of the simulated perturbers indicate a uniquely low x^, we do 
not find support for any of them. 

In order to calculate an upper mass limit above which we can 
reject perturber masses due to a high increase in;^^^, we first set 
a x^ threshold to 50.2 for TrES-1 and 28.5 for TrES-2. These 
thresholds correspond to a 90 % FAP against the linear fit us- 
ing the MPFTEST-routine. We then increased the masses for each 
semi-major axis until the;^'^ residuals reached the previously es- 
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Fig. 7. Masses of third bodies with the lowest (soHd black 
line), based on the best fit of the observed O-C values to the 
synthetic O-C values for each semi-major axis, and upper mass 
limits resulting in a FAP of 90% (dashed-dotted black line) with 
respect to the linear fit for each semi-major axis. Red lines show 
the lowest values (dashed red line) corresponding to the best- 
fit masses and the value of the linear fit is indicated by a 
horizontal dashed-dotted red line. The upper plot is for TrES- 
1 and the lower plot for TrES-2. Mean motion resonances are 
indicated by vertical dotted lines. 



tablished threshold, which gave us the upper mass limit. Masses 
above the threshold can be rejected with a FAP higher than 90 % 
(Fig. 13. 



3.3. Search for a sinusoidal transit timing variabiiity 

In order to search for any sinusoidal periodicity in the data with 
periods on the order of the campaign duration and shorter, we 
fitted amplitude and phase of a sinusoidal function to the O-C 
residuals after subtracting the linear fit with the parameters of 
TableE] A sinusoidal transit timing variation might be caused by 
an exomoon around the transiting planet ( feartoretti & Schneider 
119991: iDovle & Deeg l 120041: ISimon et aUBoOTt iKippingi 120091) . 
For this fit we stepped through dififerent trial periods between 
1500 days (6.6 x lO^^ cycl/d) and 0. 1 days (10 cycl/d) for TrES- 
1 and for TrES-2 between 600 days (1.6 x 10"^ cycl/d) and 0.1 
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Fig. 9. Results of the sine fitting as function of trial frequen- 
cies. The red shaded areas show possible frequencies which 
might have been introduced by a moon. Top: TrES-1, bottom: 
TrES-2. 



days (10 cycl/d) in steps of 0.00001, 0.0001, 0.001, 0.01 and 0.1 
cycl/d. For each trial period we then logged the x^ residuals at 
the best fitting amplitude and phase. 

From Fig.|9] showing the result of the;i^^ sine fit, we see that we 
lack outstanding peaks of low x^ values for TrES-1 and TrES- 
2, being it very unlikely that a real sinusoidal signal has more 
than one period. However, we also calculated the FAFs of the 
best sinusoidal fit against the linear polynomial, as described in 
Sect. l3.2l and show the results in Table|4]for comparison. We note 
that the best sinusoidal fit for TrES-1 with a period of 16.7 days 
and an amphtude of 25 s (Fig. (TO] upper graph) does not give 
a significant improvement against the linear polynomial, having 
a lower x^ residual but a high FAP, but for TrES-2 we found a 
good sinusoidal fit with a FAP of 1 % for a period of 0.2 days 
and an amplitude of 57 s (Fig.[TOl lower graph). 

4. Discussion 

For TrES-1 we obtained eight new O-C values and made use of 
10 previously published values of which we removed two ob- 
servations with a high probability of containing star spots. The 
standard deviation of the O-C values is ~33 s an d the maxi- 
mum O-C deviation from the established ephemeris (IWinn et al.l 
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Table 4. Comparison of different residuals for TrES-1 (upper values) and TrES-2 (lower values). 
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12007 ") is -60 s. Similarly, for TrES-2 we observed five transits 
and made use of th ree additional published mid-transit times 
jHolman et al.ll2007l) . The transit timing estimation at epoch 
for TrES-2 had an error greater than 60 s and was hence re- 
moved. After correcting the ephemeris, the standard deviation 
of this O-C diagram is ~43.9 s and the maximum O-C deviation, 
is 66.5 s. 

In order to search for a light time effect, we fitted different poly- 
nomials. Assuming the validity of the quadratic function for 
TrES-1, we can use the quadratic term to estimate a possible 
acceleration and its corresponding 1 -cr error of^ = -1.1 + 1. Ox 
10"''* s"', the negative value indicating a decelerating system. 
Comparing this deceleration to the acceleration for the solar sys- 
tem of a few 10 ''^ s ' (Zakamska & Tre maindllOOSl) and CM 
Draconis of ~10 s ' (Deeg et al. 2008i) , the acceleration of 
TrES-1 could be several order of magnitudes larger , but is not 
well c onstrained given the current data. Following iDeeg et akl 
(|2008|) . the minimum mass of a third body causing a given ac- 
celeration in dependence of its lateral distance is given by 

^> 438.26 (6) 

We note that this is independent of the mass of the acce lerating 
system . Using the distance to TrES-1 of 157 pc ( So zzetti et alj 
|2004|) and above value for ay, we obtain the minimum mass of 
a possible third body as a function of angular separation from 
TrES-1 of 




This relation allows the identification of nearby objects found 
in any future high-resolution imaging; that is, any possibly 
found object has to fulfill Eq. [T] in order to be a potential 
source of the observed acceleration. For example, for a third 
object at a distance of 0.05" the object's mass has to be at least 



0.09 Mq in order to cause the previously mentioned acceleration 
(^ = -1.1 + 1.0 X 10-'4). The angular separation of 0.05" 
translates into an orbit for a third object of 7.8 AU, correspond- 
ing to a period of 21.8 years. This stellar object would have 
caused a semi-amplitude in the radial velocity measurements 
of 1.0 km s However, the observation span of the available 
radial velocity measureme nts is only 49 days dAlonso et al.l 
120041; lLaughhnetaLll2005b . which is 0.6 % of the exemplary 
third object's period. Therefore, this object would not have 
been detected in existing radial velocity measurements. This 
example also shows that the obtained value for the acceleration 
is reasonable. 

A second approach has been assuming perturbing objects at 
nearby orbits. Then for TrES-1, perturbers with a masses >1 
at the 2:3 and 2:1 MMRs are very unlikely due to the high;t'^ 
residual with a high FAP against the linear fit. However, none 
of the simulated perturbers indicated a uniquely low peak. 
Therefore, we do not find support for any perturber, but we 
established upper mass limits above which we would obtain a 
FAP of 90 % against the linear fit. For TrES-2, we obtained 
some low peaks; none of them in a low-order MMR. For 
the best-fit simulated O-C we obtained a perturber of 18 M® 
and 0.049 AU. Due to the residuals below the one from the 
linear fit, none of the perturbers could be excluded and again 
we established upper mass limits for perturbers in the TrES-2 
system. 

ISteffen & Agoll(l2005l ) did an analysis of transit times of TrES-1, 
but they found no convincing evidence for a second planet. 
We also found general l y a h igher x'' n^ar to the low-order 
MMRs. .Steffen & Agoll (|2005|) gave upper limits for additional 
planets in the system TrES-1, whereas we give additionally 
the best-fitting mass for any orbital distance. In general, their 
upper mass limit for additional planets decreases closer to the 
transiting planet, similar to our findings in Fig.|7] We also found 
that additional planets with masses above 10 M® at the MMRs 
are very unlikely, which is consistent with the results from 
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Fig. 10. Best fit sine function (solid line) and observed O-C val- 
ues with error bars for TrES-1 and TrES-2, phased with period 
of that fit. This period is 0.06 cycles/day or 16.7 days for TrES-1 
and 4.8 cycles/day or 0.21 days for TrES-2. 



ISteffen & Agoll (l2005h . indicating that at the low-order MMRs 
and near zero eccentricity, the mass of the additional planet has 
to be below 10 M®. 

Assuming sinusoidal transit timing variations, we note that we 
find a good O-C fit for TrES-2 with a FAP of 1.1 %. However, 
since we find several good periods beyond the best one of 0.2 
days, we conclude that we need more observations and maybe 
with higher precision in order to confirm any one of them as a 
possible exomoon. Whereas for TrES-1, we found no evidence 
that a sinusoidal function improved the fit significantly with 
respect to the linear on e. 

Sartoretti & Schneider' (Il999l) . iDovle & Dee^ (|2004|) . 
Simon et al. (2007) and recently [Kipping (2009) discussed 
the possibility to detect moons around extrasolar planets 
using a timing offset induced by the wobble of the planet 
around the planet-moon barycenter. For a given timing offset 
6t the mass o f the possible moon M„, can be estimated using 
dSartoretti & Schneider.,1999,) : 



napMp 
umPp 



6t, 



(8) 



where ap is the semi-major axis of the planet, Mp the mass of 
the transiting planet, qm semi-major axis of the moon and Pp the 



period of the planet. A moon around the transiting planet should 
have a semi-major axis which is between the Roche limit VRoche 
and the Hill radius rnm (,Kipping.2009,) . The Roche hmit is given 
by: 



Roche - Rp\^ — 
Pm 



1/3 



(9) 



where pp is the density of the planet and p,„ the density of the 
moon. As can be seen from Eq.|9] if p,„ > 2pp, the Roche limit is 
inside the planet. This is most likely the case for a rocky moon 
around a gas giant planet. Therefore, we do not further consider 
the Roche limit as lower orbit limit, but the planetary radius. On 
the other hand, the Hill radius is given by: 



rHiii = ap 



Mp 
3mI 



1/3 



(10) 



where M, is the mass of the star. Using the radius of the transit- 
ing planets, i. e. RpjrES-i - 1.081 Ry or RpjrES-2 - 1-272 Ry 
and Eq.[TO] we obtain possible orbital periods for a moon around 
the transiting planets of 4.2 h < P,„ < 15 h for TrES-1 and 3.9 h 
< P,„ < 12 h for TrES-2. If we consider only the corresponding 
period ranges in Fig. |9] (red shaded areas), we also find no clear 
peak there for TrES-1, but for TrES-2 the best peak of 0.2 days 
is in the range for an exomoon. Using Eq.|8]we obtain a possible 
moon mass o f 52 for the best - fit arn plitude of this peak. 
We note that iBarnes & O'BrienI (|2002|) gave an analytical ex- 
pression for an upper mass limit of possible moo ns. Generally, 
the m ore massive the moon the shorter its lifetime. [Lubow et al.l 
119971) estimated that close-in orbiting planets spin down 
into synchrono us rotation very quickly. Therefore we use the 
procedure from lBarnes & O'BrienI (l2002h to estimate the upper 
mass limit of a moon around TrES-1 and TrES-2. We obtain 
an upper mass limit around TrES-1 on the order of ~10 *' 
and around TrES-2 on the ord er of ~10~^ Mm, s i milar to one 
obtained for HD 209458b by iBarnes & O'BrienI (l2002l) . This 
means that moons with masses greater than 10~* will not 
have survived until now. This upper mass limit is clearly way 
below our detection limit for masses causing transit timing 
variations in the system TrES-1 and TrES-2. The closest moon 
at this upper mass limit would cause a timing amplitude of the 
order of 10"*' s, which is not detectable even from space by 
several orders of magnitude. 

Tra nsit dura tion investigations on TrES-2 performed by 
Mis lis & Schm itt (2009) indicated a decrease of 3 min. between 
2006 and 2008 in the duration. However, in our transit timing 
analysis, we do not find non-linear deviations of O-C times 
with a similar magnitude. While we can not provide any firm 
detection, we can put upper mass limits which are consistent 
with our observations for TrES-1 and TrES-2, respectively. 
However, due to the gaped data, it might have been possible 
that we missed important points in the observed O-C diagram, 
like e. g. a transit timing measurement at the highest amplitude 
at the mean-motion resonances. We need at least about 3 years 
of continuous transit observations and with several transit 
observations per year in order to avoid the missing of windows 
of high-amplitude O-C deviations from perturbing bodies, with 
the maximum amplitudes being observable only during a few 
transit events, see Fig.|6] But even for low-mass perturbers such 
amplitudes should be easily measurable from ground at mean- 
motion resonances. On the other hand, outside the mean-motion 
resonances, we need higher precisions than obtainable from 
ground to detect the small amplitudes caused by a perturbing 
body of planetary mass, something that may be expected from 
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the current satellite missions in operation. Any transit timing 
variation in TrES-2 may be expected to be confirmed in the near 
future by observations done with the Kepler Space Telescope. 
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